(N 

o 

(N 



On 



> 

00 

oo 



% 



E.M. Ovsiyuk* 

Quasi-plane waves for a particle with spin 1/2 on 

the background of Lobachevsky geometry: 

simulating of a special medium 



-gj February 27, 2013 

■*— > 

S . 

Abstract 

In the paper complete systems of exact solutions for Dirac and Weyl 
equations in the Lobachevsky space H3 are constructed on the base 
of the method of separation of the variables in quasi-cartesian coordi- 
nates. An extended helicity operator is introduced. It is shown that so- 
^o ■ lution constructed when translating to the limit of vanishing curvature 

t^J" , coincide with common plane wave solutions on Minkowski space going 

in opposite z-directions. It is shown the problem posed in Lobachevsky 
f~^ , space simulates a situation in the flat space for a quantum-mechanical 

C^l ' particle of spin 1/2 in a 2-dimensional potential barrier smoothly rising 

to infinity on the right. 



It is known that in the field theory of elementary particles, the basis 
of plane wave is of the most use. However, in presence of a curvature, 
any common plane wave solutions do not exist. Therefore, of a special 
interest are examples non-Euclidean spaces in which some analogues of such 
solutions can be constructed. In the paper [I], it was shown that in the 
Lobachevsky space there are such solutions for particles with spin 0; also 
see the books by Gelfand-Graev-Vilenkin [3] , [5] . An analog of plane waves 
in a space of constant positive curvature was studied by Volobuev |2j. The 
later treatment of this problem was given in [3] . Solutions of the plane wave 
type for Maxwell's equations have been considered in [6]-[9]. In |10j . the 
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problem of constructing solutions of the Dirac equation in the Lobachevsky 
space was studied on the base of the method of squaring; in particular, it 
was pointed out the possibility of constructing solutions of the Dirac plane 
wave studied starting with Shapiro's scalar waves . In this paper we will 
construct a complete basis of solutions of the plane wave type for Dirac and 
Weyl particles in the Lobachevsky space, applying the method of separation 
of the variables in a special system of quasi-cartesian coordinates closely 
related to horospherical coordinates. 

To understand the physical meaning of the system under consideration, it 
should be mentioned that Lobachewsky geometry simulates a medium with 
special constitutive relations. The situation is specified in quasi-cartesian 
coordinates (x,y,z) was treated in [9]. Exact solutions of the Maxwell equa- 
tions in complex 3- vector form, extended to curved space models within the 
tetrad formalism, have been found in Lobachevsky space. The problem re- 
duces to a second order differential equation which can be associated with 
an 1-dimensional Schrodinger problem for a particle in external potential 
field U(z) = Uoe 2z . In quantum mechanics, curved geometry acts as an ef- 
fective potential barrier with reflection coefficient R = 1; in electrodynamic 
context results similar to quantum-mechanical ones arise: the Lobachevsky 
geometry simulates a medium that effectively acts as an ideal mirror. Pen- 
etration of the electromagnetic field into the effective medium, depends on 
the parameters of an electromagnetic wave, frequency ui, kf + k 2 , and the 
curvature radius p. 

In the present paper, that analysis will be extended to the case of parti- 
cles with spin 1/2, described by equations of Dirac and Weyl. The general- 
ized spinor plane waves can find application in the analysis of the behavior 
of fermions particles on cosmological scales, or in simulating special media 
affecting the spinor particles. 

1 On the solutions of the Schrodinger equation 

In the Lobachevsky space-time parameterized by quasi-cartesian coor- 
dinates 

dS 2 = dt 2 - e~ 2z {dx 2 + dy 2 ) - dz 2 ; 

the element of volume is given by 

dV = y/—g dxdydz = e~ z dxdydz , x,y,z G (— oo, +oo) . 

The magnitude and sign of the z are substantial, in particular when referring 
to the probabilistic interpretation of the wave functions 

dW =| ^ | 2 dV =| tf | 2 e~ 2z dxdydz . 



Let us describe some details of the parametrization of the space by coor- 
dinates (x,y,z). It is known that this model can be identified with a branch 
of hyperboloid in 4-dimension flat space 



2 2 2 2 2 / o o 

u -u 1 -u 2 -u 3 = p , u = +VP + u • 

Coordinates x,y,z are referred to u a by relations 

u\ = xe~ z , u 2 = ye~ z , 

u 3 = ±[(e z -e- z ) + (x 2 + y 2 )e- z }, 

u Q = l -[(e z + e- z ) + {x 2 + y 2 )e' z \. (1.1a) 

It is convenient to employ 3-dimensional Poincare realization for Lobachevsky 
space as an inside part of 3-sphere: 

ft = — = 1 5 ./ 5 . o ftft < +1 • (1-1&) 

Quasi-cartesian coordinates (x, y, z) are referred to qi as follows 

2x 2y 



91 ~ o | .,9 I .2, I 1 > 12 



x 2 + y 2 + e 2z + 1 " x 2 + y 2 + e 2z + 1 

x 2 + y 2 + e 2z - 1 

93 = 2 , 2 i 2s i 1 • L1C 

z + y + e zz + 1 



Inverses to (1.1c) relations are 



ft 92 2 \/l-g 2 ,, n n 



1-^3 1-93 1-93 

In particular, note that on the axis 9i = 0, 92 = 0, 9 G (— 1, +1) relations 
(l.ld) assume the form 



n -■ ,1 + 93 
x = 0, y = 0, e = x , 

1-93 



that is 



93 — 


->+l, 


e z - 


->■ +CX) , 


2 — 


->■ +00 ; 


93 - 


->-l, 


e z - 


^+0, 


2 — 


-7- — 00 . 



(1.2) 



Schrodinger equation in Riemannian space |llj 



\hdt9 



1 

2M 



i h 



(-=d ky /=jj + eA k )(-g* l )(ihd l + eA i ) 

v-5 



* 



in quasi-cartesian coordinates (1.1a) takes the form 

d 2 ,. a 2 



•* d T 



2M 



e 2z — + e 2z — + e 2z — e~ 2z — 
dx 2 dy 2 dz dz 



$. 



The variables are separated by the substitution \& = e lEt / n e lk i x e tk ^v f(z): 

f2 



^-2^- + e-e 2z (k 2 + k 2 ) 
dz 2 dz 2J 



f(z) = 



(1.3a) 



where a dimensionless quantity used e = 2MEp 2 /fi , p - curvature radius 
of the space. Elementary substitution / = e z ip(z) in equation (1.3a) gives a 
Schrodinger-like equation 



dz 2 



+ e-l-(k 2 + k 2 )e 2z <p(z 







with potential function 



2\2z 



U(z) = l + (kf + ki)e 



[1.3b) 



[1.3c) 



Note that the probabilistic interpretation of the wave function after the 
transformation to <p> reads 



dW =| * | 2 dV =| if | 2 dxdydz . 



(1.4) 



An easily interpretable physical solution for e > 1 is the following: on 
the left we have the superposition of two waves, falling from the left and 
reflected. On the right behind the barrier, the wave function must sharply 
decrease to zero. 

It should be noted that the case k\ = 0, &2 = is special: the equation 
(1.3a) is very much changed because the potential function disappears 



3?- 2 S + «)/M = ° 



/ 



3 (l±i N /<PT)z 



<P 



,{±i-Je-L)z 



(1.5) 



and the function <p is a solutions of the type of ordinary plane wave. 



Let us turn to the general case and in eq. (1.3a) introduce the variable 



'k 2 + k%e z = Z , Ze(0,+oo); 
the equation takes the form 



£ l d 



dZ 2 ZdZ z 2 



+ « 1)/(Z) = 0; 



(1.6) 



with the help of a substitution / = y Z F, one can remove the term with 
the first derivative 



d 2 e - 3/4 



+ 



1 F{Z) = . 



dZ 2 Z 2 

This form makes it easy to find the asymptotical behavior of solutions 
(z — > — oo) Z — > , 



(z — > +oo) Z — > +oo , 



ip ~ e 



±i\/e— 1 z 



F ~ e 



±z 



/ = ^e ±z 



ip r^ e z ' exp 



2 c 2 



±v/ fc i + A;^ e 



(1.7) 

We now turn to the construction of exact solutions of (1.6) in the entire 
range of coordinate z. We seek solutions in the form of f(Z) = Z e F(Z); 
equation (1.6) gives 

d 2 F dF 



+ \{B 2 - 1)Z -B(1-2A) + 



A(A-2) + e 



F = 0. 



(1- 



At A, B chosen according (for definiteness, we take the minus sign before 
the root in the expression for A; assuming e > 1) 



A= 1-iy/e-l 
the equation (1.8) is simplified 



B 2 = 1 



d 2 F dF 

z -J^ + ( 2A ~ l + 2BZ ) -J^-B(1-2A)F = 0. 
dZ A dZ 



(1.9) 



;i.io) 



In (1.10), let us make another change Z = y/2: 

d 2 F dF \ 

with B = —1 it is an equation for the confluent hyper geometric function 



y — + (2A-l + By) — + B(A--)F = 0, (1.11) 



d 2 Y , ^dY 



c = 2a , a = A - 1/2 = 1/2 - We- 1 , 

f{Z) = y a+l l 2 e-yl 2 Y{y). (1.12) 

We use two pairs of linearly independent solutions [12] 

Yi = *(o, 2a, y) , y 2 = y^^l - a, 2 - 2a, y) 

and 

y 5 = *(a,2a,y), Y 7 = e^(a,2a, -y) , (1.13) 

These pairs of solutions are related by Kummer linear relations [12] 
_ r(l - 2a) T(2q - 1) 

y5 "r(i-a) yi+ r(a) y2 ' 

T(l-2a) r(2o-l) 

^7 = -^ f li - \, , ' Y 2 , 1.14a 

1 (1 — a) 1 (a) 

which after multiplication by y a + l l 2 e~ y ' 2 take the form 

_ r(l-2o) r(2a-l) 

75 r(i - a) Jl r(o) /2 ' 

f _ r(l-2a) r(2a-l) 

Note that the solutions Y\ and Yi describe the case of negative z — > — oo 
wave with the asymptotic behavior 

z -+ -oo , (y — >0) 

h ~ V a+1/2 = (V A;? + h*\ 1 l ^ e z e- l ^ lz , 



a+1/2 ,.l-2o = 9 / t 2 , 1,2 ^ p +iA/i=T ^ 



/ 2 ~ V a+l/z V l - za = 2Jk 2 + k 2 e z e w * . (1.15) 



Thus, for example, the function Y5 (and the related ip§) at negative 
z — > —00 behaves as a superposition of two plane waves according to 



¥?5 



+ 



T(l - 2a) 

" r(i-a) 
r(2o-i 



2Jk\ + kl 



\-iJ7^I 



2Jk\ + kl 



1+iWe-i 



e -iVe-lz + 



+iy/e— 1 z 



(1.16) 



We define the reflection coefficient as the square modulus of the ampli- 
tude ratio in a superposition of plane waves 



M_ e - iV7 ^ z + M+ e +ivCT z 



R 



M_ 



M 4 



i? 



r(l - 2a) r(a) 



We take into account 

1 - 2a = +2We - 1 
a = 1/2-iVe-l , 



r(2a - 1) T(l - o) 

2a - 1 = -2zVe + l 



(1.17a) 



l-o= 1/2 + iVe-l, 



then 



R 



T(+2i^F=T) 



T(-2i^/F^l) 



r(i/2 - iVe^I) 



1 . 



(1.176) 



r(i/2 + iv^^^) 

We find the behavior of I5 at large y. Using the known asymptotic 
relation |12j 

Y 5 = V{a,c,y)~ y - a , 



we get 



~ [2Jk( + kie 



+00 , h = y a+1 ' 2 e- y / 2 Y 5 ~ y 1 /2 e -j//2 
1/2 

-* ' '•-<■"- > -- -0. (1.18) 



2 „z 



exp 



rCi ~\~ rZc\ C 



exp 



Thus, the solution /s describes the expected situation: wave going from 
the left is reflected with probability 1 on the effective barrier; behind the 
barrier the solutions sharply decrease to zero. It is easy to find the critical 
point, after which wave function must sharply decrease 



2\ „2z 



e - 1 = (jfef + ki) e 



zq = In , 



1 



rC-i ~t~ /Co 



1.19a) 



in the usual units, this critical point is described by the relation 



<2mEp 2 /h 2 



1 



The solution ft for y 



26 -""v W+^V 

+00 goes to infinity 



(1.196) 



ft ~ y a+l/2 e -y/2 e y y -a 



y +m e +y/2 



However, it is of no clear physical interpretation. 

However, it is easy to interpret such a solution ft. Indeed, far on the 
right behind the barrier, taking into account the asymptotic formula [12] 



we get 



+00 , 



+00 



HAC,y) = ^ e yy A - c . 



T(A) 



h 



y i/2 e y/2 



+oo L 



1.20) 



that is, for z — > +00 is a real function with infinite probability density. Far 
left is a plane wave propagating to the left. 

Next we consider the analogue of this situation for particles with spin 
1/2, described by the relativistic Dirac equation, when analysis is much more 
complicated. 



2 The Dirac equation in the space H 3 , separation of variables 

We start with the general covariant form of the Dirac equation [11] 

d , 



i l a ie? 



d 1,1 



-0 e (a). 



"(a) Q x a ' 2 K ,/=gdx c 

In the coordinate system (1.1) we use the tetrad 



/?? 



*(x) = 0. (2.1) 



» 



1 














e z 














e z 














1 



eq. (2.1) takes the form 



ff (7 (7 1 Ci 

^d-t + ^ Z d- X + ^ Z 0- y + ^[d-z 



1 I + im 



* = 0. 



(2.2) 



(2.3) 



Note that the addition of —1 about the operator d z can be removed by 
substituting ^f = e Zr ip. The following three operators idt, id x , id y commute 
with the wave operator: solutions can be searched in the form 



q,e,ki,k 2 _ „-iet J,k\x e ik 2 y 



h{z) 

Hz) 

h{z) 



(2.4) 



Using the Dirac matrices in spinor basis, from (2.3) we find e equations for 
fi(z) 

-ie f 3 - ik x e z U ~ k 2 e z fi - ( j- - 1) h + im fi = , 

-ie f A - ik ± e z f 3 + k 2 e z f 3 + ( — - 1) / 4 + im f 2 = , 
-ie h + ihe z h + k 2 e z f 2 + (jL-l)f 1 + i m f 3 = 0, 



-ie f 2 + ik t e z f\ - k 2 e z f t - ( — - 1) f 2 + im / 4 = 



(2.5) 



There is a generalized helicity operator which commutes with the oper- 
ator of the wave equation: 



S--(rVvV + e* 7 3 7 1 !- + 7 1 7 2 (7--l) 



d_ 

Ox 



dy 



d_ 

" dz 



(2.6) 



Using the substitution (2.4) in the eigenvalues equation £ ^ = p *$> we 
obtain 

he z f 2 - ik 2 e z f 2 - i( — -l)f 1 = pf\ , 
he z h + ik 2 e z f 1 + i(-^--l)f 2 = pf 2 , 
he z U ~ ik 2 e z f 4 - i{ -r- - 1) h = Ph , 

he z h + ik 2 e z f 3 + i(-^--l)f 4 = pf 4 . (2.7) 

dz 

From equations (2.7) and (2.5), considered together, it follows a linear ho- 
mogeneous system with respect to fi 



-ie f 3 - ipf 3 + im f x = 



-ie f4 - ip f4 + im f 2 = , 

-ie /i + ip /i + im f 3 = , 

-ie h + ivh + imf 4 = 0. (2.8) 

We we find two values for the p and the corresponding restrictions on the 
functions ff. 



±V ^T^ ; h = LJP. hi U = L-Ef 2 . (2.9) 



Thus, we have three continuous quantum number e, k±, k 2 and one dis- 
crete, which distinguishes the values p = ±Ve 2 — Tn 2 . In view of (2.9), from 
four equations (2.5) we arrive at two equations for /i, f 2 

(^-l-ip)fi + e z (ik 1 + k 2 )f2 =0, 

(^-l + ip)f2-e z (ik 1 -k 2 )f 1 =0. (2.10) 

Note the symmetry of the equations with respect to change 

/i=>/ 2 , P=>-P. (2-11) 

It is convenient to obtain solutions of similar equations in the flat space 

( ^ - ip) /i + ( ifci + to) h = , 



so that 



{^- z +ip)f2-{iki-k 2 )h =0, (2.12) 

^ + e 2 -m 2 -k 2 1 -ki)f 1 = 0. (2.13) 



Here, there exist two independent solutions (let k% = + \/e 2 — m? — k\ — k 2 ) 

f (i) _ p +ik 3 z ,(i) _ (+ifc 3 ~ JP) +ifc32 

/x " C ' h ~ ik x + k 2 6 

ff) = e -«s* , /W = _ (~ jfe 3 ~ «P) e -^ . (2 . 14) 

ifci + fc 2 

The sign before k% determines the direction of the wave propagation, the 
sign of p defines the state of polarization. Generalized analogue of these 
solution is to be investigated in the hyperbolic space H%. 

10 



3 A special case of the waves along the z-axis 

There exists a special case when k\ = 0, k 2 = 0: 



yW(t,z) =e 



-id 



h{z) 
h{z) 

Hz) 

U(z) 



The equations change substantially (see (2.5)) 



d 



■it fa ~ ( t ~ 1) h + im h = 

dz 



(3.1) 



-it h + ( j- z ~ 1) h + im h = 
-*e /i + ( ^ - 1) /i + im / 3 = 

-*e / 2 - ( ^- - 1) / 2 + »m / 4 = 

diagonalization of the operator X gives (see (2.7)) 

(|-i)A-*A, (|-DA- 



-ip/2 



(^ - 1) /3 = ip/3 , (^ - 1) /4 = -ip/4 



(3.2) 



(3.3) 



Considering equations (3.2) and (3.3) together, we arrive at the linear system 
(e + p) /3 - m /i = , (e + p) f± - m f 2 = , 



(e - p) h ~ m / 3 = 
From this it follows 

( e 2_ ;p 2 ) 2_ m 4 = 

m 



(e - p) /2 - m U = 



(3.4) 



p = ±v e 2 — m 2 



f - e " P f 

74 — 72 



/2 



m e + p 

Equation (3.2) reduce to 

( £ - 1 - ip) h = 

az 



f - e_;p f 
73 — /l 



m 



in 



e + p 



f 



i • 



(3.5) 



h = Ci eV 



-ipz' 



11 



( Tz ~ 1 + ip) h = ° => 

Solutions more simple to interpret are 
d = 1 , C 2 = , tf e00 P(i, z) 



h = C 2 e z e~ 



ipz 



(3.6) 



d = , C 2 = 1 , * e00p (t, . 



1 


e-p 

m 




1 



e-p 
m 



e z e +i P z 



(3.7a) 



e z e-~ ipz . 



(3.76) 



Obviously, the factor e z in the solutions will be compensated when con- 
sidering any bilinear structure of the wave functions (with their subsequent 
multiplication by y/—g dxdydz). 

4 Construction of solutions in the general case 

Let us turn to (2.10) and introduce a new variable Z: 



kf + k^e z = Z, Ze(0,+oo) 



d 



z s- i -*j* + *v£3£ a =°- 



4- i+ *A + V^t A =°- 



It is convenient to define a parameter 



k 2 + ifci 
^2 — ik\ 



then 



Z—-l-ip)f 1 + Ze^ a f 2 = 0, 



z dz~ 1 + ip i f2 + Ze ~ mfl = ° • 



(4.1a) 



Note that the factor e ia can be removed by changing the notation for f 2 



e +ia f 2 



h\ 
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then the system takes the form 

z-^ - 1 - ip) fi + Zf 2 = 0, 

Z-^-l + ip)f 2 + Zf 1 = 0. (4.1b) 

From (4.1a) we get two second order differential equations for f\ and fi'- 

Note the symmetry between the equations: they are transformed into each 
other when changing p — > —p. Furthermore, it should be noted that unlike 
the case of flat space (see (2.13)), here second-order equations for the func- 
tions /i and /2 depend explicitly on the first degree of p, that is depends on 
the state of polarization of spinor waves. 

Considering eq. (4.2a), let us use a substitution f\(Z) = Z e BZ F\(Z): 



+ 



,2 1W , no ,A ^ , (A-ip-l){A + ip-2) 



(B l - l)Z + 2B(A-l) + 



Fi = . (4.3) 



Z 
At A and B chosen according 

A = +ip + l, -ip + 2, B = ±l, (4.4) 

(4.3) becomes simpler 

Z 4t^ + (2A-2 + 2BZ) ^- + 2B(A- l)F 1 = . (4.5) 

dZ z dZ 

The resulting equation with another change of Z = y/2 will transform into 

y-j^ + (2A-2 + By)-£±+B{A-l)F 1 =0. (4.6a) 

When B = — 1, it coincides with an equation for the confluent hypergeomet- 
ric function with parameters (for definiteness let it be A = +ip + 1) 

d 2 $ d$ 

y — -jj- + (c - y) — a<J> = 0, a = +ip , c = 2a = +2ip. (4.66) 

ay ay 

13 



Two linearly independent solutions are |12j 

Fl 1 \y) = $(a,c,y), 

F[ 2) (y) = y x - c $(a - c + 1, 2 - c, y) . (4.6c) 

Consider the equation (4.26). Using the above-noted symmetry, we ob- 
tain 

h = y a ' +1 e- y/2 F 2 (y) , a' = -ip , c' = 2a' = -2ip , 

F 2 (1) = d>(a',c',y), 

F® =y 1 -*$(a' -<t + l,2-<t,y) . (4.7) 

It is convenient to employ one independent parameter a: 

h = y a+1 e-y' 2 F l {y), 

F^(y) = ^(a,2a,y), 
F? ] (y) = y 1 - 2 "^ -a,2- 2a, y) ; (4.8a) 

f2 = y- a+1 e-^ 2 F 2 (y), 
F 2 {1) (y) = <S>(-a,-2a,y), 

F 2 {2) (y) = y 1+2a $(l + a, 2 + 2a, y) . (4.86) 

The functions f\, f 2 (note that before now we did not find possible nu- 
merical factors at them) must be related by the first-order operators (see 
(4.1)) 

4- 1 -°) /l -f e+ta/2 = ' 

»5-i + «)A-|«-A-0i 

these relationships can be translated to the functions F\, F 2 , which results 



m 



^±-~(F 1 + y~ 2a e +ta F 2 ) = , (4.9a) 

ay 2 

^--\{F2 + y +2a e-* a F 1 )=0. (4.96) 
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These equations relate functions in the following pairs 

F^ (y) if } (y) , F[ 2) (y) if } (y) . 

For each pair one should find the relative coefficient of the two functions. 
Let us substitute expressions for F^ (y) and F 2 (y) with some numerical 
coefficients 

F™ (y) = rj 1} $(a, 2a, y) , F 2 (2) (y) = rf y 1+2a <D(a + 1, 2 + 2a, y) 

in eq. (4.9a). Performing the necessary differentiation 

rj 1} [$ (a + 1, 2a + 1, y) - $ (a, 2a, y)] - rf ] e +ia y $ (a + 1, 2a + 2, y) = 0, 

and transformation using the relations for contiguous functions, we find the 
relative factor 

r [V = 2 e +ia (2a + 1) rf ) . (4.10) 

The same result can be obtained by substituting the expression for F[ (y) 

(2) 

and F 2 (y) in the equation (4.96). The resulting ratio can also be obtained 
by using the expansions of solutions near y = 

*i (1) = r«(i + ^ + ^j/ + -, 4%) = r? v 1+2a ; 



the equation 



^^l F (l)_l y -2a e+iajF (2) 



dy 2 * 2 
gives 

(1) / 1 a + 1 1 1 \ +ia (2) n 

ri U + 2oTT y "2-2 y J- e r * y = ° => 

r?— J__-c+ fa rJ ,) =0. 
1 2(2a + l) 2 

Now let us substitute expressions for FJ (y) and -F 2 (y) of (4.8) 

F 1 (2) (y) = rf ) y 1 - 2a $(l-a,2-2a,y), 

F 2 (1) (y) = 4 1} (y)$(-a,-2a,y) (4.11a) 

in the equation (4.96). Performing the differentiation 

r[ 2) [$ (1 - a, 1 - 2a, y) - $ (-a, -2a, y)] - 
15 



-r^ 1} e' ia y$(l-a,2-2a,y) = 

and the corresponding transformations, we arrive at 

r {2) =2e- ia {l-2a)r i l ) . (4.116) 

Thus, the construction of two linearly independent solutions of (4.1): 
/l=r ( 1 ) e ^/V +1 $(a,2a, 2/ ), 

h = rf e~y/ 2 y 2+a Ha + 1,2 + 2a, y) , 
rj 1 ) = 2 e +ia (1 + 2a) r {2) ; (4.12a) 

/2 = r ( 1 ) e ^/V a+1 $(-a,-2a,y), 
/l=r f) e ^/y- a $(l-a,2-2a,y), 

r f) = 2 e~ ia (1 - 2a) r£ 1} . (4.126) 

We simplify the formulas by setting r 2 = 1 , r 2 = 1; the result is 
/ f 1 = M + e-y/ 2 y 1+a ®(a,2a,y), 

f 2 = e-y/ 2 y 2+a Ha + l,2 + 2a,y), 
M + = [2e+ ia (l + 2a)] ; (4.13a) 

// h = M_ e~yl 2 y 2 - a $(l - a, 2 - 2a, y) , 

f2 = e^ 2 y 1 - a H-a,-2a,y), 

M_ = [2e~ ia (l-2a)] . (4.136) 



Remind that a = ip = ±i\/e 2 — m 2 ; the sign of p is associated with the 
polarization state of the spinor waves; types / and II are supposed to be 
associated with the directions of wave propagation: to the left or to the 
right. 

Let us consider asymptotic properties of the solutions. The coordinate 
y = 2y/k'f + k\ e z tends to (see (1.2)) 



1, y = 2Jk( + k 2 2 e z — ► 0, 



+1 , y = 2^Jk 2 + k%e z — ► +oo. (4.13c) 
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The behavior of solutions in boundary point y — > (z — > — oo) 
I 

1+ip 



2+ip 



l+a _ i\j / o. u2 , j,2 c £ 



h = V 2+a = Uy/% + k% e 2 ) W ; (4.14a) 



1—ip 



II 



f2 = y 1 - a =[2yjk 2 1 +k 2 2 e z \ (4.146) 

For large values of y one should use the asymptotic formula [12] 

y^+oo, ^A,C,y) = J ^ e yy A - c . 

So we get (for z — > +oo, y — > +oo) 

T(a) T(a + 1) 



(4.15a) 



r(l — a) F(— a) 

(4.156) 

To conclude this section we consider the limiting process in the con- 
structed solutions (4.13a), (4.136) to the case of the flat space. This will 
allow a better understanding of the obtained results in the Lobachevsky 
space. 

To this end, we first need to go to the usual dimensional quantities: 

z 3 McR ER 

z= r> m = — ' e = ^T' 



P 



+Ve 2 - m 2 = +R JE 2 /c 2 h 2 - M 2 c 2 l% 2 = Rp , 
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a = ip = iRpo, c = 2a = i2Rp . 



y = 2sjkl + k 2 e z = 2RK ± (1 + ^ + ...) — > 

Let us consider the solutions (4.13a) 

f 2 = e -y/ 2 y 2+a Ha+ 1,2 + 2a,y) 
M+= [2e +ia (l + 2a)] ; 
taking into account 



2i?K, 



(4.16) 



a 1 
-c y= 2 y 



RK\ 



1 o(a + l) 2 1 l/2( 1 / 2 + 1 / c ) 2 

-y = ^r — /-, rm — y 



2!c(c+l) tf 2! (1 + 1/c) 



2! (jR ^ 



la(o + l)(o + 2) 2 _ 1 1/2(1/2 + l/c)(l/2 + 2/c) 2 



3! c(c + l)(c + 2)" 2! (l + l/c)(l + 2/c) 

1 



we get 



,-y/2 



-w/2 



e 

-i?A" 



-itK"i 



3! ( ^ )3 



$(a,2a,y) 
$(a + l,2a + 2,y) 



(4.17) 



,R/fi 



R^i 



and further 



/ /i 



M+ (2iLK"_Le 



2\l+iRp .. ^3P0 



/2 



(4.18) 



=► {2RK ± e z ) 2+iRpo ~ e" 3P0 . 

Similarly, we find 

// / 2 = e -^ 2 y 1 -«$(-a, -2a, y) ~ e -»«» , 

/i = M_ e^V^^l - a, 2 - 2a, y) ~ e~ ia;3Po . 

(4.19) 

We may conclude that solutions of the type I (in curved model ^3) 
provide us with extension for the flat waves in Minkowski space of the type 
g+jfcz. w h ereas solutions of the type II represent extension for the flat waves 
in Minkowski space of the type e~ . 
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5 The case of the Weyl neutrino 



Let us restrict the above analysis to the case of 2-component Weyl neu- 
trino. It is enough to specify the Dirac equation in the spinor basis (see 

mi) 

i (J a {x) [ d a + E a (x) ] £(x) = m r]{x) , 

i a a (x) [ d a + S a (x) ] r}(x) = m £(x) 

and impose relation m = 0, thus we will obtain the equations of the Weyl 
neutrino (wave function rj(x)) and antineutrinos (with the wave function 

Substitution for rj is 

h(z) 



V 



e,fci,fc 2 



e ~iet ^ikxx e ik 2 y 



h{z) 



(5.1) 



Note that neutrino wave function has only three quantum numbers e, k\, k 2 . 
The system of equations after separation of variables has the form 



d 
-ie h ~ ihe z fa - k 2 e z / 4 - (— - 1) / 3 ■ 

-ie U - ihe z f 3 + k 2 e z f 3 + (^ - 1) / 4 



0. 




(5.2) 



or (replacing / 3 on hi and / 4 on h 2 



1 + ie) hi + e z (iki + k 2 ) h 2 = 



1 — ie) h 2 — e z (iki - k 2 ) hi = 

az 



(5.3) 



One can use the solution obtained above (for the system (2.10)), replacing 
everywhere p on — e. Thus, we get two linearly independent solutions of the 
system (5.3) 



where 



hi = r[ 1} e~ y/2 y a+1 $(a,2a,y) , 




h 2 = rf e-^y +a $(a + 1,2 + 2a, y) , 




r (l) =2e +m (1 + 2a)r (2) . 


(5.4a) 


/l2=r ( 1 ) e ^/V a+1 1>(-a,-2a, 2 /), 




fo = r{ 2) e^V -0 *^ - a, 2 - 2a, y) , 




r [ 2 ) = 2e- ia (l-2a)4 1) , 


(5.46) 


a = — ie , c = 2a = — 2ie . 


(5.4c) 
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6 On the non-relativistic Pauli approximation 

We will carry out a procedure of the non-relativistic approximation di- 
rectly in the separated equations for relativistic Dirac case. To this end, we 
turn to equations (2.5) 

-ie f 2 + ihe z h - k 2 e z A - ( — - 1) f 2 + im / 4 = , 

dz 

-ie U - ihe z / 3 + k 2 e z f 3 + ( — - 1) f 4 + im f 2 = , 

-ie h ~ *he z / 4 - k 2 e z f 4 - ( -j- - 1) / 3 + im /i = , 

dz 

-ie h + ihe z f 2 + k 2 e z f 2 + ( -£■ - 1) h + im f 3 = . 



We introduce new functions 

/i + /3 _ r /i — /3 

2 ~ J ' 2i ' 



5 , 



(6.1) 



^±A = F , !l_l± = G . (6 .2) 

From (6.1), one obtains equations for /, F, g, G : 

(-^-l)G-e z (ih -k 2 )g+(e-m)F = 0, 
{^--l)F-e z (ih - k 2 ) f - (e + m) G = , 

{ Tz~ l)9 + eZ {ih + fc 2) G - (e - m) / = , 

(-£■ - 1) / + e z (ifc! + fc 2 ) F + (e + m) g = . (6.3) 

dz 

Now we should make a formal change e =^> m + E (this is equivalent 
to separating the rest energy by a factor e~ rmt ); as a result we get 

(4- - 1) G - e 2 (ifo -fc 2 ) 5 + EF = 0, 
az 

(— - 1) F - e z (ih - k 2 ) f - (E + 2m) G = , 
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,d_ 

v dz 



( l) g + e z (ik l + k 2 )G-Ef = 0, 

dz 

1) / + e z (ih + k 2 ) F + (E + 2m)g = 0. 



(6.4) 



The condition for the applicability of the nonrelativistic approximation 
is the following relation E + 2m ~ 2m; which results in 

(4- ~ 1) G - e z (ih -k 2 )g + EF = 0, 
dz 



G= — 

2m 



' 1) F - e z (ih - k 2 ) f 

dz 



' dz 



l)g + e z (ik 1 + k 2 )G-Ef = 0, 
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1 
2m 



(^--l)f + e z (ih + k 2 ) F 
dz 



(6.5) 



Excluding two small components of g, G, we arrive at two equations for 
the big components / and F: 



^ 1 -2^ + l-e 2z (kl + k 2 2 ) + 2mE 
dz 2 dz \ i zj 

^-2^- + l-e 2z (k 2 + k 2 ) + 2mE 



F-e z (ih-k 2 )f = 0. 
f + e z (ih + k 2 )F = 0. 



We remind that (see (2.9)) 

fi(z) + f 3 (z) 1 + A 



f 

F 

A- 



2 2 

f2(z) + h(z) 1 + A 



h{z) , 

Hz) , 



€± P i /~2 2 

, p = ±\/ e z — m* 



in 



±y/(E + m) 2 - m 2 « ±V2mE 



(6.6) 



(6.7a) 
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sign ± correspond to two different polarizations of the non-relativistic elec- 
tron. Therefore, equation (6.6) can be represented as follows (remember 
that now the parameter p defined by the nonrelativistic expression (6.7a)): 



^-2^- + l-e 2z (k 2 + k 2 ) + 2mE 

^-2^- + l-e 2z (k 2 + k 2 ) + 2mE 
dz 2 dz \ l z) 



f 2 -e z (ik 1 -k 2 )f 1 = 0, 



fi + e z (ik 1 + k 2 )f 2 = 0, 



(6.76) 



These equations can be written as (for definiteness, let p = +\/2ME) 

f 2 -e z (ik 1 -k 2 )fi = 0, 



{ Tz ~ l ~ ip)i Tz " l + ip) " &2z (fc? + k2) 



{^- z -l-ip){^- z -l + ip)-e 2z {k 2 + k 2 ) 



h + e z (ik 1 +k 2 )f 2 = 0. 



We take into account (2.10) 



d 



( j- z ~ l ~ *P) h + e "( ik i + k ?)h = , 



(^-l + ip)h-e z (ik 1 -k 2 )f 1 =0. 



Thus we arrive at equations 



d 



(^ - I - ip) (e z (ih - k 2 ) h)- 



-e z (ih - k 2 ) (— - 1 - ip) h - e z (ih -k 2 )h=0. 



dz 



1 + ip) (-e z (ih + k 2 ) f 2 )+ 



d 



+e z (ih + k 2 ){—-\ + ip) f 2 + e z (ih + k 2 )f 2 = 0, 



(6.9) 



It is easy to see that these two identities of the form = 0. 

Thus, the solutions found above for the relativistic Dirac equation in the 
approximation p ~ \/2mE are solutions of two-component Pauli equation. 
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7 On representation of solutions in terms of Bessel functions 

Returning to the basic system (4.1) in the form 

Z— -1-ip) \Jk 2 - ifafi + Z\Jk 2 + %k\h = , 

Z— -1 + ip) \Jk 2 + ikif 2 + Z\Jk 2 - ikifi = , 
let us express it in terms of new functions 

s/k 2 -ihfi = e z ipi , \Jk 2 + ihf 2 = e z ip 2 ; (7.1) 

as a result we obtain 

Z dZ ~ ip ) Vl + ZiP2 = ° ' 

Z lZ + ip ) W + Zip^O. (7.2) 



Let us translate them to the new variable x = iZ = i\Jk\ + k 2 e z : 



d \ 

x- ip I ipi — ixip 2 = , 

x— + ip ) (p 2 - ixipi = . (7.3) 

From (7.3) it follow two second-order equations 

d 2 p^ + ip 

(-r-2 +i + — —m = o, 

dx A x z 

< * + 1 + 1^1 )V2 = . (7.4) 

dx z x A 

Separating the factor y/x : (f\ = y/xF±, ip 2 = y/xF 2 , we arrive at the two 
Bessel equations 

\ dx z x dx x z J 

f d 2 Id p 2 - ip- 1/4 \ 

\ dx 2 x dx x 2 J 
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Since we are (especially) interested in the case of the Weyl neutrino 
(with helicity —1), we will continue to consider in detail the case of negative 
p = —\/e 2 — m 2 . To avoid new notation, we will make a minus sign in front 
of p in (7.5), so that 

fd 2 | 1 d | x | pJ»-fe-i/4 \ 

V dx 2 x dx 



p 2 — ip - 


-1/4 


X 2 




p 2 + ip - 


-1/4 



^ + ^ + 1 + '-^^!'- ft 



where p = +y/e 2 — m 2 > 0. 

Thus, for polarization states with p = +\/e 2 — m 2 > 0, the functions 
Fi(y),F2(y) satisfy the Bessel equations [13] 



# + I< +1 _«£±|0>!Wo, 

ax z x ax x z J 

v = -ip- 1/2 , F\ = J+v(x), J- V (x) . (7.7a) 

\ax z x ax x z J 

H = -ip + 1/2 = v + 1 , F 2 = J +M (x), J_ M (x) . (7.76) 



Transition to polarization states with p' = —y/e 2 — m 2 > is achieved 
by replacement in (7.7) positive p to negative p' = —p, which results in 



z/ = +%p - 1/2 = -/j, , // =+ip + l/2 = i/ + l = -v. (7.7c) 

Since we have to follow explicit form of both functions Fi,F 2 (up to a 
relative factor) , let us return to the first order equations (now in the variable 
x = iZ); with notation v = —ip — 1/2 it reads 

x— -v ] Fi = -ixF 2 , (7.8a) 

dx J 

x— + u + l)F 2 = -ixF 1 . (7.8b) 

dx J 

From the above it is known that functions F\(x),F 2 (x) satisfy the Bessel 
equations (7.6). Let us recall the well-known recurrence formulas for the 
solutions of Bessel equation |13j - write them in a convenient form 

x- v J F v (x) = -xF v+ i(x) , (7.9a) 
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x— - iA F_ v {x) = +xF„ u _ l {x) ; (7.96) 

here by F± v can be understood either of the Bessel functions: J± u , or Hankel 
functions H± u , H± UJ or Neumann functions N± v {x). 

Comparing eq. (7.9a) with eq. of (7.8a), we find two types of solutions: 

in Bessel's functions 

J F{(x) = J+„(x) , F 2 (x) = -i J +{u+1) (x) ; 



II F{ T (x) = J- V {x) , F 2 H {x) = +i J- (v+ i)(x) 



in Hankel's functions 



/ F{{x) = H\ v (x) , Fi(x) = -iHl (u+1) (x); 



(7.10) 



II F/'(x) = H% v {x) , Fi\x) = -i H 2 +{v+1) {x) ■ 






(7.11a) 



I' F 1 {x) = Hl u {x), F 2 (x) = +iH 1 _ (u , l) (x) 



-U\^J } A W) — T' -"_(V +1 )\ 



II' F x {x) = H 2 _ v (x) , F 2 (x) = +i H 2 , +1) (x) ; 



(7.116) 

note that H^_ v (x) = e lv *H~(x), so the primed cases I', II' coincide respec- 
tively with II, / and by this reason will not be considered below. 
And in Neumann functions 

/ F((x)=N +v (x), F 2 I {x) = -iN +{u+1) {x)- 
II F{ I (x) = N_ u (x), F 2 H {x) = +iN_ {v+1) {x). (7.12) 



First, let us detail solutions in Bessel's functions [13]. 
In the region z — > — oo, x — >• iO , 



F((x) = J +U (x) = (-)" = (iXY e-^e-*/ 2 , 
Fi(x) = -i J +(u+1) (x) = -i (i\y +1 e-^e^l 2 
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II 



Fi T (x) = +i J_ ( „ +1) (x) = +i {i\)- v - 1 e + ^e- z ' 2 ; 



(7.13) 



here and below we use the notation 



VW+M 



In the region z — > +00, x — > +ooi, using the knows asymptotic formula [13] 

T ^ H^ ( t 1 ^\ 
J v [x) ~ \ / — cos \x — (u + -)— , 



Tr.r 



we get 



2 / pn 

J +U {z ->■ 00) ~ V ~~^ cos * (^ + -y 



2'2 



1 



2vriX 



e+^/ 2 e x 



J- U (z — )• OO) 



2 . ./ y P7T\ 

— — sin i A 

ivrX V 2 / 



Thus, the solutions (7.10) behave 

2 — >■ +00, x = iX — > +ioo, 



1 



2-iriX 



e -pn/2 e x 



(7.14) 



I F{(x) = J +V (x) ~ J^ e^ e x 



F 2 7 (z) = -i J +(l/+1) (s) ~ +a/^ e +p7r/2 e* 



2vriA" 



// F( x {x) = J- V {x)~i 



1 



2vriX 



e - P -/2 e x 



F 2 // (x) = +iJ_ (l/+1) (x)~-y^ 



2mA" 



e -^/2 e X 



(7.15) 



Let us consider solutions (7.11a) in Hankel's functions [13]. They are 
determined in terms of J± v (x) as follows 



SVO.VK 



H l v {x) = +— z (e-^J +v {x) - J- V {x)) , 
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Hl{x) = -^— (e +i ™J + „(x) - J_„(x)) . (7.16) 

From here it is easy to set behavior of the Hankel functions for small x: 

z — > — oo, x — > iO , 

HUx) ~ +^— (e- iu1r (i\Ye- ipz e- z / 2 - {i\y v e ipz e +z ' 2 
sinz/7r V 



i 



sinra 



~ iun (iX) v e~ ipz e- z/2 



Hl{x) ~ — (e +iu7T (i\) u e~ ipz e- z/2 - (i\y v e ipz e +z/2 ) 

sinz/7r V / 

e +iu *{i\y e- ipz e- z ' 2 , 



smvir 

i 



smvir 

(7.17a) 



Hl +1 (x) ~ + -^ — (e- i{v+1 ^{i\y +1 e- ill,z e +z/2 - {i\y v ' l e +ipz e - z/2 

smiv + 1 )7r V 



sin(f + l)-7r 

(iA)"^ 1 e +ipz e 



' ■■■ -v-l +ipz -z/2 



Hl +1 (x) ~ 



sin(z/ + l)7r 



sin(z/ + l)7r 



' {i\T v ~ l e +ipz e~ z / 2 



sin(z/ + 1)tt 

(7.17b) 

Therefore the solution of the system in Hankel functions behave as fol- 
lows 

z — > — oo, x — > iO , 

I F((x) = HlJx) — e -i™{i\y e - ipz e- z l 2 , 

sin vk 

Fi(x) = -iH\ {v+1) {x) ~ ~* (a)— i e +fr« e -«/* ; 

^ v T ; sm(i/ + 1)7T 



// F/ 7 (x) = H 2 {x) e +iu7r (i\ye- ipz e~ z / 2 , 

sin vk 
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F/ J (x) = -iHl (v+1) (x) ~ -_i__(iA)- 
TlT) sin(f + IJ7T 



y-1 +ipz -z/2 



'e ' ^ e 



Behavior of the Hankel functions of z — > +00 is given by 



(7.18) 



H l( x ) ~ V — ex P 

7TX 




•I T 1 "/ In 

+1 \x--(y + -) 




H u {x) ~ a / — exp 
V vrx 

so that 

z — » +00, x = iX — > +ooi 



-1 (x--(i/ + -) 



(7.19) 



Hl(x) 



mX 
H 2 v {x) 



exp 






ivrX 



e -W2 e -x 



mtA" 



exp 



- UX --(, + - 



ivrX 



e +P*/2 e +X 



^H-lM 



iirX 



exp 



+/ -(ix-|(i/ + i + i; 



~ — 1 



mX 



e -pn/2 -x 



Hi+M 



iirX 



exp 



.j[ix-~(v + i + h 



+1 



mX 



e + P 7r/2 e +X 



(7.20) 



/ F/(x) = H\ v {x) 



mX 



e -W2 e "* 



F 2 ! (x) = -zff' 1} (x 



fvrX 



e -P*/2 e -x 



// F^(x) = ^(x) 



ivrX 



e + ^/ 2 e+ x 
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F»{x) = -i H*(x) ~ */ — e + ^ 2 e +x . (7.21) 



2 
l +(H-l)W - y -^ 

Finally, let us consider solutions in terms of the Neumann functions. 
These functions are as follows expressed in terms of Bessel functions [13] 



COSVK J v (x) - J- U (x) 

N u (x) = 



saxvn: 



= UX)-<XX™J-U(X) 

sin vk 
In the region z —> +oo, (x —> +ooz) we use the asymptotic formulas 



v„U,^/-^sin(*X-(, + I)f 



p ~X-ia _ p +X+ia _• 

sin(iX - a) = -e x e +ia 

y ' 2 2 



we get 



N +U (z -> oo) ~ +J -!—<?« I* e x 



JV +1 , +1 (z -► oo) V ^ eW2 e * 



A_„<: -, x)- \/^J e ~ W2 e X 



#_„_!(* -+ oo) ~ y ^ e -W2 e x (7>23) 

Thus, the solution (7.12) in the functions of the Neumann behave at 
z — > +oo as follows: 



/ Fl{x) = N +V {x) ~ i\ — — e 



1 



C W2 P X 



2-kzX 



e 



Fi(x) = -i N +{u+1) (x) ~ +i yilI e W2 e x . 



// /-;"(,)- v_ |/ ,,.,. M /_L_ e -f»r/a e * 



F/ J (x) = +» M ( „ +1) (s) ~ - x /-l^ e -^ 2 e x . (7.24) 
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Now, using the definition (7.22), we find the behavior of the Neumann 
functions at z — > — oo. 

cosz^7r J„(x) - J- V (x) 

N v {x) = : 

smz/7r 

cosz/vr {iXf e^e" 2 / 2 - (iX)~ u e +ipz e +z / 2 







sin^7r 










rsj 


COSZA7T 

sinz/7T 


- (i\) u e 


-jpz g -; 


s/2 

5 






AT , <^\ - 


cos(z^ 


+ l)vr J„ 


+i0) 


- J_ 


-i/- 


-i(x) 



sin(f + l)-7r 

cos(^ + l)vr (iA)^ 1 e - ipz e +z / 2 - {iX)" 1 " 1 e +i P z e~ z / 2 
sin(z/ + l)-7r 

' (iX)"^ 1 e +ipz e~ z / 2 , 



sin(z/ + l)7r 

J^(x) - COSZ/7T J_^(x) 
N_ l/ (x) = ; 

smz/7r 

(iA)" e-^e" 2 / 2 - cosz^vr {iX)~ v e +i P z e +z ' 2 
sinz/7r 



smz/7r 



AU-i(x) 



— (iA/2)" e^e"*/ 2 , 

IVTT 

Jv+i{x) - cos(u + l)ir J_ v _i(x) 



sin(z/ + l)7r 

(iAf +1 e -^ z e +z / 2 - cos(z/ + 1)tt (iA)-^ 1 e+^e-*/ 2 
sin(i^ + l)-7r 

_ -CQS(Z/+1)7T ( - A) -,-l e +i P z e -z/2 _ 



sin(z/ + lW 



Consequently, the solution (7.12) in the functions of the system behave 
Neumann z — > — oo as follows: z — > — oo, x — >■ iO , 



/ F/(x) = tf + „(s) ~ ^L (iA)" e-^e-/ 2 , 



itf(s) = -i i\r +( „ +1) (s) ~ -i \ (iX)-"- 1 e +^ e -/ 2 ; 

v ' Sin(i/ + 1)7T 
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// F^ix) = NJx) (iXY e~ ipz e~ z l 2 , 

sin vn 

E?(x) = +i N_ (v+1) (x) ~ +* ~ C 7 ( " + 1)7r (a)— 1 e+^e-/ 2 . 
v ; sm(i/ + IJ7T 

(7.25) 

The results can be collected into a table. 

a = —p 

Bessel Hankel Neumann 

z — > — oo z — > +00 z — > —00 2: — > +00 z — > —00 2 — > +00 



e -*p^ e +jY 



Z[ e~ ipz e +x Z[ e~ ipz e~ x Z[ 

Z{ e +x Z{ e +ipz e~ x Z{ e +ipz e +x 

Z II q e +X Z II e -ipz e +X Z II e -ipz e +X 

Z II e +i[z e +X Z II e +ipz e +X Z II e +ipz e +X 

The most interesting are solutions of the type / in Hankel functions: 

z -> -00 , Z{ ~ e" ip2 , 4 ~ e +ipz , 

z -> +00 , Z[ ~ e" x -*• , Z| ~ e" x -> , 

which means that the problem posed in Lobachevsky space can simulate a 
situation in the flat space for a quantum-mechanical particle of spin 1/2 
(Dirac electron and Weyl neutrino) in a 2-dimensional potential barrier 
smoothly rising to infinity on the right. Electromagnetic field behaves itself 
similarly - see in 



Conclusions and acknowledgement 

In the paper complete systems of exact solutions for Dirac and Weyl equa- 
tions in the Lobachevsky space H3 are constructed on the base of the method 
of separation of the variables in quasi-cartezian coordinates. An extended 
helicity operator is intoduced. It is shown that solution constructed when 
translating to the limit of vanishing curvature coincide with common plane 
wave solutions on Minkowski space going in opposite z-directions. Electro- 
magnetic field behaves itself similarly. 

Author is grateful to V.M. Red'kov for encouragement and advices. 
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